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ABSTRACT 

A  method  based  on  a  theory  of  experimental  Optimal  Design  (OD)  was  developed  in  Russia 
to  select  the  best  combination  of  wavelengths  so  an  optically  remote  spectral  sensor  may 
optimally  estimate  oceanic  chlorophyll  concentration.  It  gives  the  number  of  spectral  bands,  their 
center  wavelength,  and  each  band’s  spectral  width  to  make  the  best  estimate  of  chlorophyll 
concentration  in  a  fixed  observation  time.  A  hyperspectral  sensor  is  convenient  because  it 
eliminates  the  need  for  variable  optical  filters.  An  ideal  sensor  would  have  an  infinite  number  of 
infinitesimally  narrow  spectral  bands.  Optimal  designs  are  given  for  both  ideal  and  real  sensors. 

The  designs  were  computed  using  ocean  radiance  spectra  simulated  by  a  Monte  Carlo  model  for  a 
range  of  chlorophyll  concentrations.  The  computed  OD  was  tested  for  robustness  over  a  wide 
range  of  experimental  conditions.  The  methods  of  color  index  and  principal  component  analysis 
were  also  applied.  The  optimal  design  method  gives  more  accurate  chlorophyll  estimates. 

1.0  INTRODUCTION 

An  important  task  for  the  modem  remote  sensing  oceanographer  is  the  accurate  estimate  of  in  situ  chlorophyll 
(or  other  in  situ  parameters)  from  spectral  measurements  of  water  leaving  radiance.  A  method  was  developed 
(Kozlov,  et  al.,  1992,  Levin,  et  ai,  1994,  1996,  and  1997,  and  Zolotukhin,  et  aL,  1997)  to  optimally  select  spectral 
channels  for  the  remote  determination  of  an  in  situ  ocean  parameter.  This  method  was  based  on  a  theory  of 
Experimental  Design  originally  suggested  by  V.  Kozlov,  (1974  and  1983).  It  has  a  few  advantages  over  the  present 
methods  of  color  index  and  principal  component  analysis  used  for  chlorophyll  estimation.  The  wavelengths 
considered  in  these  latter  methods  are  often  qualitatively  selected  and,  thus,  do  not  guarantee  an  optimal  use  of 
spectral  channels.  Spectral  bands  generally  have  the  same  width,  e.g.,  5,  10,  or  20  nm,  with  the  result  that  optimal 
spectral  width  is  not  treated.  Additionally,  noise  found  in  real  sensor  data  is  not  addressed.  In  contrast,  the  Optimal 
Design  (OD)  algorithm  searches  for  the  ^y  optimal  spectral  channels,  giving  center  wavelengths  and  band  widths 
while  accoimting  for  thermal  and  shot  noise  (Levin,  et  ai,  1994,  1997,  and  Kozlov,  et  a/., 1989). 

The  OD  algorithm  uses  as  input  the  first  two  moments  of  the  joint  distribution  of  remotely  measured  ocean 
radiance  spectra  u(A)  and  the  parameter  of  interest  ©  (here,  chlorophyll  concentration).  The  spectra  u{X)  are 
expressed  as  the  number  of  primary  electrons  in  each  of  the  sensor’s  spectral  bands  per  unit  time.  As  noted,  the 
parameter  0  may  be  the  concentration  of  phytoplankton,  sediment  or  some  other  optically  active  matter.  The 
moments  of  the  distribution  (u(X),  0)  are:  <0>  and  cr^,  the  mean  and  variance  of  =  <u(A)>,  the  mean  value 
of  observed  spectral  photoelectrons  u(X)\  K{A^X)  =  cov  (w(>l),  w(A')),  the  auto-covariance  function  of  w(A);  and 
q(X)  =  cov  (u(XX  0),  the  covariance  between  u(X)  and  0. 

Given  the  statistics  on  (w(>l),  0),  an  optimal  design  §  may  be  computed  which  determines  the  number  of 
independent  channels  r  of  the  sensor,  the  spectral  transmittance  function  x,(>l)  for  each  channel,  and  the  time 
periods  allotted  for  measurements  in  each  channel.  The  f/  are  constrained  by  Zr,  =  /,  i.e.,  the  overall  observation 
time  is  fixed.  Since  x,{l)  is  the  spectral  transmittance  fimction  of  the  /-th  channel,  then  0  <  Xi(A)  <  1.  Each  channel 
may  consist  of  multiple  spectral  bands.  The  signal  observed  by  the  optimally  selected  filter  jc,(1)  is  represented  as 
the  number  of  primary  electrons  in  the  /-th  channel  in  the  time  period  r,. 
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Kozlov,  et  aL,  (1989)  showed  that  the  best  linear  estimate  0  of  the  parameter  0  is  given  by 

©  =<0>+/D"'r'r.  (1) 

with  the  residual  variance  of  this  estimate,  considered  to  be  a  measure  of  design  optimality  or  “goodness,”  given  as 

h{^  =  a^-q^D-^q,  (2) 

where  q  =  (q,,...qj,  q,  =  \qiX)x,iX)dX,  D  =  K  +dmg(J,  +v)7’-'.  K  KiJ  =  \\K(,X,X’)x,{X)Xj{X)dXdX' , 

fi  =  T  =  diag(/^),  Y  =  CViv^r)^*  the  noise,  independent  of  the  signal  level.  They  also 

proved  that  the  optimal  design  consists  of  no  more  than  two  spectral  channels  (r  =  2),  taking  on  values  of  only  0  or 
1,  i.e.,  there  will  be  only  two  filters,  Xx{X)  and  X2(A),  each  consisting  of  a  number  of  rectangular  transmission  bands. 

The  algorithm  was  originally  used  by  Levin  (1996,  1997)  and  Zolotuchin  (1997)  to  estimate  chlorophyll 
concentration  C  in  the  ocean  for  a  sensor  with  a  photomultiplier  and  a  spectral  filter  device  which  allocated  the 
measurement  time  t  between  two  spectral  channels  xi(A)  and  jc2(A)  into  periods  t\  and  t2  respectively.  The 
spectrometer  SFS-1  (Afonin,  etaL,  1987),  with  a  rotating  filter  set,  is  such  a  sensor;  sequentially  sampling  different 
spectral  bands. 

In  this  paper,  we  apply  the  OD  theory  to  the  problem  of  ocean  remote  sensing  with  a  hyperspectral  sensor 
system,  which  eliminates  the  need  for  variable  optical  filters.  Accordingly,  noise  is  considered  that  of  a  semi¬ 
conductor  detector,  and  radiation  is  received  simultaneously  by  both  channels  over  the  entire  measurement  time, 
corresponding  to  the  condition  0.5/,  where  /  is  the  total  measmement  time  used  in  the  OD  calculation,  not 

the  sensor  frame  time  (which  is  actually  /i  or  The  spectral  width  of  the  hyperspectral  sensor,  although  small,  is 

finite,  thus  imposing  a  limit  on  the  realization  of  an  “ideal”  optimal  design. 

2.0  SIMULATION  OF  INPUT  RADIANCE  SPECTRA 

The  OD  method  requires  the  use  of  a  large  database  of  upwelled  radiance  spectra  related  to  in  situ 
chlorophyll.  The  collection  of  such  a  database  is  prohibitively  expensive.  Thus  the  approach  taken  was  to  simulate 
the  spectra  for  the  database  using  a  Monte  Carlo  oceanic  radiance  model.  Atmospheric  light  and  scattering  were 
made  negligible  by  only  modeling  the  data  that  would  be  obtained  from  shipboard  or  low-flying  airborne 
observations. 

The  OD  method  uses  the  joint  distribution  statistic  (i/(A),  Q  for  spectral  photoelectrons  and  chlorophyll.  The 
spectrum  w(A),  the  spectral  density  of  the  primary  electrons  per  unit  time  (e.g.,  [p  *  s‘*]),  is  obtained  using  the 
radiance  spectral  density  U  (A).  Lt  (A)  is  the  sum  of  spectral  densities  of  the  water-leaving  radiance  just  above  the 
surface,  L  ^A),  and  of  the  radiance  reflected  from  the  sea  surface,  L^A),  with  units  of  [W  m'^  sr"'  p*)].  Thus 


u(X)  =  L,  (X)  A  t(X)  CO  q{X)e^  (X), 

(3) 

L,(X)  =  LXX)  +  LXX), 

(3a) 

whereof  [m^],  i(A),  w  [sr],  and  //(A)  are  the  sensor  parameters,  i.e.,  aperture  (entrance  pupil)  area,  transmittance, 
receiver  solid  angle,  and  photo-detector  quantum  efficiency  (QE).  e{X)  is  the  photon  energy  {e  [J]  =  2x10* V  A  [p]). 
L  'is  related  to  the  reflectivity  coefficient  p  =  TdJEd  (L  is  the  spectral  density  of  the  water-leaving  radiance  just 
beneath  the  sea  surface)  as 


LXX)  =  Cs7f^p(X)E(AX 


(4) 


where  E{X)  [W  m'^  |i  *]  is  the  spectral  density  of  solar  irradiance  of  the  sea  surface  and  Cs  =  0.533  (Austin,  1974) 
takes  into  account  the  radiance  change  through  the  water-air  interface.  The  sea  surface  radiance  for  a  rough 
surface  is  obtained  from  p/A)  =  Substituting  for  L  and  is  given  as 

LiX)  =  7f^E{X)  [0.533p(>l)  +  psim  .  (5) 

We  consider  only  a  cloudless  sky.  The  irradiance  E{A)  is  a  sum  of  irradiances  from  direct  sunlight  Ef'XX)  and 
from  dififiisive  sky  light,  A),  assumed  to  be  uniform.  For  each  irradiance,  water  leaving  and  surface  reflected, 
the  corresponding  coefficients  take  the  different  values  p/O  and  pt^\  In  the  general  case, 

p  -  aff'""  +  (1  -  ps  =  apt"  +  (1 "  d)pt^,  (6) 

where  a  is  the  ratio  of  direct  sunlight  to  the  total  sea  surface  irradiance. 

The  coefficients  pt"  w^ere  computed  for  nadir  viewing  direction  by  Mullamaa's  formula  based  on  the  Cox- 
Munk  surface  slope  distribution  (Mullamaa,  1964).  They  are  functions  of  wind  velocity  v,  azimuth  angle  and 
solar  zenith  angle  6|).  Table  1  shows  several  examples  of  calculated  values  of  pt"^.  For  a  uniform  sky  in  nadir 
direction,  pt^  ^  0.02  is  virtually  independent  of  wind  velocity  (Mullamaa,  1964). 


Table  1.  Surface  Radiance  Coefficients 


The  values  of  surface  radiance  coefficients  p/'  (averaged  over  azimuth  <p^ 
for  direct  sunlight  in  the  nadir  viewing  direction 

Solar  zenith 
Angle  6>° 

Wind  velocity  v  (m/s) 

2  5  10  15  20 

0 

0.42 

0.20 

0.10 

0.072 

0.054 

20 

0,036 

0.064 

0.061 

0.051 

0.043 

45 

0.0001 

0.0016 

0.008 

0.014 

0.018 

The  coefficients  and  //^for  nadir  viewing,  accounting  for  particle  and  molecular  scattering,  were 
calculated  (Levin,  1997)  as: 


dir  ^  l-?7^  + 0.8(1 +  0.8//o^);7,  \ 

2(l+//o)  a  +  b^  ’ 

=  (0.27  +  0.01  t]J)  bt  /(a  +  bt),  (8) 

where  a  and  bi,  are  the  absorption  and  backscattering  coefficients,  //6  =  cos[arcsin(sin^«)],  t)„  =  biJbt,  and  bi^  is 
the  backscattering  coefficient  of  optically  pure  water.  a(X)  and  bb{X)  were  obtained  as  functions  of  concentrations 
of  chlorophyll  (Q,  yellow  substance  (f)  and  sediment  particles  (X)  through  the  use  of  the  models  of  Prieur  and 
Sathyendranath  (1981)  and  Sathyendi^th,  et  al.  (1989).  Following  Sathyendranath,  et  al.  (1989),  the  ranges  of 
C,  A" and  y  were  selected  to  be  C  =  0.01  -  100  mgm■^A"=0.01  -  10  m”',  and  7  =  0.01  -  1  m"'. 

One  thousand  water  samples  of  varying  C,  X,  and  Y  were  chosen  for  Monte  Carlo  modeling  by  randomly 
sampling  the  above-mentioned  C,  X,  and  Y  distributions.  The  logarithm  of  each  concentration  was  assumed  to  be 
uniformly  distributed  over  the  range  given  above.  The  solar  angle  ^  was  distributed  uniformly  within  the  limits  of 
40-60°  (at  smaller  solar  angles  the  amount  of  sun  glint  in  the  receiver  view  increased  sharply).  The  wind  velocity  v 
followed  a  Rayleigh  distribution  with  mean  velocity  <v>  =  10  m/s,  and  the  azimuthal  angle  ^  was  uniform  within 
the  limits  of  0-180°.  The  solar  irradiance  density  E(X)  for  various  0o  was  taken  from  tables  (Shifirin,  et  al.,  1959), 
and  the  ratio  a(6b,  A)  from  experimental  data  (Austin,  1979). 


For  each  combination  X,  Y,  C,  do,  v  and  the  spectra  were  computed  using  Eqs.  5-8  and  the  lOP  model 
of  Prieur  and  Sathyendranath  (1981)  and  Sathyendranath,  et  al.,  (1989).  A  database  containing  1000  simulated 
radiance  spectra  L^A)  with  a  step  of  5  nm  was  obtained. 

3.0  THE  SENSOR  PARAMETERS 

Since  the  OD  is  calculated  from  the  photoelectron  counts  in  the  sensor  detector  elements,  sensor  parameters, 
including  noise,  that  affect  the  photoelectron  count  must  be  taken  into  consideration.  The  computations  were 
conducted  for  a  hyperspectral  sensor  with  aperture  area  A  =  5.7  lO’^m^  (lens  diameter  2.7  cm),  field-of-view 
0)=  =  2.4x10'''  (2/3=9.94°  is  a  total  viewing  angle  over  all  320  sensor  elements),  t  =  10.5  ms,  QE  =  0.6. 

A  sampling  of  the  other  sensor  parameters  required  for  the  OD  computations  that  vary  with  spectral  detector 
site  is  given  in  Table  2.  This  table  gives  the  detector  spectral  channel  number  (No),  the  center  wavelength  A,  of  that 
channel  in  run,  the  bandwidth  AA  of  the  channel  in  run,  the  total  transmittance  t  of  all  the  optics  and  detector 
coatings,  the  gain  g  of  the  channel  in  photon-to-electron  conversion  units,  and  the  signal  independent  noise  vin 
number  of  electrons.  The  latter  includes  shot  noise  of  the  thermal  background  radiation,  digitizer  noise,  electronic 
readout  noise,  and  offset  correction  error.  Parameters  presented  in  Table  2  are  typical  of  a  hyperspectral  sensor  with 
an  InSb  detector  array.  In  addition,  there  is  a  gain  noise  factor  F  =  1.3  that  is  the  same  for  all  elements. 


Table  2.  The  Sensor  Parameters 


No 

X 

AX 

r 

g 

y 

No 

X 

AX 

T 

g 

y 

1 

400.610 

3.07 

0.39 

1.95 

1919 

28 

507.781 

4.94 

0.56 

1.79 

2 

403.906 

3.11 

0.39 

1.97 

1899 

29 

512,89 

5.04 

0.55 

1.76 

1877 

3 

407,229 

3.15 

0.41 

1,97 

1905 

30 

518.119 

5,16 

0.56 

1.74 

1882 

4 

410.582 

3.20 

0.42 

1.99 

1899 

31 

523A12 

5.28 

0.56 

1.70 

1881 

5 

413.967 

3.24 

0.43 

1.99 

1901 

32 

528.953 

5.4 

0.56 

1.67 

1881 

6 

417.388 

3.28 

0.44 

2,01 

1902 

33 

534.588 

5.53 

0.56 

1.63 

1876 

7 

420.847 

3.34 

0.44 

2.01 

1900 

34 

540.323 

5.58 

0.56 

1.61 

1862 

8 

424.348 

3.39 

0.44 

2.03 

1896 

35 

546.222 

5.8 

0.57 

1.56 

1872 

25 

493.121 

4.62 

0.55 

1.87 

1882 

52 

674.06 

8.88 

0.59 

1.46 

1864 

26 

497.901 

4.73 

0.55 

1,85 

1883 

53 

683.511 

9.1 

0.59 

1.46 

1867 

27 

502.785 

4.82 

0.56 

1,81 

1880 

54 

693.204 

9.33 

0.59 

1.44 

1874 

4.0  THE  OPTIMAL  DESIGN 


The  simulated  spectra  Lt(X)  and  the  sensor  parameters  noted  above  were  substituted  into  Eq.  3  to  calculate  the 
distribution  of  spectral  photoelectrons,  u(X).  The  appropriate  statistical  moments  -  <C>,J(X),  K(A,X),  and  q(X)  - 
were  then  computed.  Then  the  method  noted  earlier  (Kozlov,  et  al.,  1992,  Levin,  et  al.,  1994, 1996,  and  1997,  and 
2tolotukhin,  et  al.,  1997)  was  used  to  compute  the  “ideal”  OD  for  a  sensor  with  infinitesimal  bands  sharing 
measiu'ement  time  between  two  filter  configurations.  All  other  sensor  parameters  (A,  oa,  t(>l),  r],  t,  g,  v,  and  F)  used 
were  those  of  the  real  sensor.  Detector  shot  noise  was  computed  by  multiplying  the  square  root  of  the  number  of 
signal  electrons  by  the  gain  g  and  the  noise  factor  F.  The  resulting  Optimal  Design  for  this  ideal  sensor  is  as 
follows:  The  first  combination  of  spectral  bands,  for  the  spectral  transmittance  function  xi(A),  was  400-403, 533- 
548, 586-614,  and  686-700  nm.  The  second  combination,  for  X2(A),  was  440-459, 548-586,  and  660-667  nm.  The 
optimal  time  of  measurement  in  the  first  channel  was  t\  =  0.46r;  in  the  second  one  t2  =  0.54t. 


A  best  estimate  of  the  chlorophyll  concentration  C  is  given  by 


logC=Ao+Ai>'i  -Ajyi- 


(9) 


where  >^1  and  >*2  arc  the  signals  measured  in  the  two  channels  xi(A)  and  X2(A),  expressed  in  number  of  primary 
electrons.  The  coefficients  Aq,  Ai  and  v42,  which  are  calculated  by  the  OD  method  through  Eq.  1,  depend  only  on  the 
product  of  the  sensor  and  measurement  parameters  Aa)t=M,  but  not  on  the  values  of  the  individual  parameters. 
The  value  of  M  used  in  the  calculations  corresponded  to  the  sensor  defined  in  section  3.0  and  gave  coefficient 
values  of^o  =  0.012,  =  3.66x10“^  andyl2  =  6.46x10"^.  For  other  sensors  with  different  values  ofM,  Eq.  1  must 

be  used  again  to  calculate  the  coefficients  Ao,A\  and>l2.  The  “goodness”  of  the  estimate,  h(^  =  0.0099,  i.e.,  the 
variance  of  estimate  log  C,  is  calculated  using  Eq,  2.  This  variance  calculation  came  from  averages  over  the  entire 
set  of  chlorophyll  cases  simulated. 

Optimal  Design  for  the  real  sensor  differs  from  that  for  the  computed  ideal  sensor.  Here  the  measurement 
times  are  equal,  t\  =  ^2,  the  sensor  bands  are  of  finite  width,  varying  with  wavelength.  Therefore,  the  real  sensor 
OD  closest  to  the  ideal  uses  the  spectral  band  weights  for  the  first  and  the  second  channels  shown  in  Tables  3a  and 
3b  respectively.  The  signals  from  the  sensor  bands  in  the  first  line  of  the  tables  are  to  be  multiplied  by  the 
corresponding  weights  in  the  second  line  of  the  tables. 


Table  3a  Weights  for  Channel  1 


Sensor  band 

1 

33 

34 

35 

42 

43 

44 

45 

54 

Signal  weight 

1.0 

0.69 

1.0 

0.87 

1.0 

1.0 

1.0 

0.50 

1.0 

Table  3b.  Weights  for  Channel  2 


Sensor  band 

13 

14 

15 

16 

17 

36 

37 

38 

39 

40 

41 

51 

1.0 

1.0 

1.0 

1.0 

0.72 

1.0 

1.0 

1.0 

1.0 

1.0 

0.64 

1.0 

For  the  real  sensor,  the  Optimal  Design  gives  the  coefficients^o=  0.012,  Ax=  3.60x10“^  .42= 6.39x10"^  with  a 
“goodness”  calculated  by  Eq.  2  of  =  0.015,  the  variance  of  the  estimate  log  C.  For  comparison,  the  method  of 
color  index  was  applied  to  the  entire  simulated  data  set.  The  specific  indexes  used  were  I\  =  Z/(443)/L(550)  and 
I2  =  L(520yL(550).  The  method  of  least  squares  was  used  to  obtain  the  goodness  parameters,  which  for  the  two 
indexes  chosen  were  0. 13 1  and  0.067  respectively.  Additionally  the  method  given  by  Sathyendranath,  et  al.^  (1989) 
was  applied  to  the  data  set.  This  method  used  the  5  wavelengths  -  410, 445,  520,  565  and  640  nm  -  with  a 
triangularly  shaped  spectral  window  of  half-width  5  nm.  The  goodness  or  variance  of  the  estimate  using  this 
method  was  0.039,  In  summary,  the  goodness  values,  i.e.,  the  variances  of  the  estimate  log  C,  were  as  follows:  for 
Optimal  Design,  /i(4)  =  0.0 15;  for  the  two  color  indexes,  =  0. 13 1  and  0.067;  and  for  the  method  of 
Sathyendranath,  et  ai,  1989,  =  0.039.  Thus  the  Optimal  Design  method  provided  the  highest  retrieval  accuracy 

of  the  methods  considered. 


5.0  TEST  OF  OPTIMAL  DESIGN  ROBUSTNESS 

The  computed  ideal  optimal  design  was  tested  for  robustness  and  stability  using  more  stringent  conditions  as 
follows:  1)  The  receiver  related  characteristic M-Acot  was  taken  to  be  100  times  smaller,  and  the  detector  noise 
used  was  100  times  greater  than  in  basic  model.  2)  The  solar  angle  range  was  extended  to  ft  =  30°  to  70°  rather 
than  the  range  ft  =  40°  to  60°  in  the  basic  model.  This  meant  the  working  day  was  lengthened  with  observations 
made  closer  to  dawn  and  dusk.  3)  A  narrower  range  of  chlorophyll,  sediment,  and  yellow  substance  concentrations 
were  considered,  i.e.,  C  =  0.02  to  20  mg  m■^  X  =  0.02  to  3  m"\  T  =  0,01  to  1  m"\  Additionally,  per  Kopelevich,  et 
aL,  (1994)  a  relation  between  C  and  Y  was  introduced  where  Y  -  0.427  (-0.008  +  0.261C)  when  C  >  0.4.  And 
finally,  4)  The  Inherent  Optical  Property  model  (lOP)  of  Kopelevich,  et  ai,  (1994)  replaced  the  model  of 
Sathyendranath,  et  ai,  (1989)  and  IMeur,  et  ai,  (1981). 

Table  4  summarizes  the  resulting  impact  on  chlorophyll  estimation  under  the  just  mentioned  parameter 
changes.  The  second  column  names  the  changed  parameter.  Each  time  a  change  was  run,  the  resulting  filter 
combinations  and  weights  were  slightly  changed  from  the  base  case.  The  third  column  gives  the  new  variance  in 


the  log  chlorophyll  estimate  with  the  new  Optimal  Design.  The  fourth  column  shows  what  variance  would  be  if  the 
filter  combination  and  weights  computed  for  the  base  sensor  case  were  left  unchanged.  Note  that  for  this  case. 
Model  number  0,  the  two  variances  are  equal.  Robustness  is  demonstrated  by  noting  that  even  for  extreme  changes 
of  input  parameters  different  from  the  base  case,  the  worst  variance  in  the  fourth  column  still  only  deviates  by  a 
factor  of  1.8  from  the  new  Optimal  design.  Even  for  this  worst  case  the  goodness  or  variance  shown  is  still  from  4 
to  14  times  better  than  that  of  the  color  index  and  principal  component  methods  previously  discussed. 


Table  4.  Variance  of  Estimate 


Model 

Number 

Parameters  varied  for  robustness  test 

Variance  of  log  CMorophyll  estimate 
Changed  OD  Baseline  OD 

0 

Basic  model 

0.0099 

0.0099 

1 

Sensor  resource  M 

0.0114 

0.0150 

2 

Range  of  solar  angle  do 

0.0110 

0.0155 

3 

Range  of  concentrations  X,Y,C 

0.0096 

0.0121 

4 

lOP  model 

0.0103 

0.0188 

6.0  CONCLUSIONS 

For  data  gathered  with  remote  optical  spectral  sensors,  the  optimal  design  algorithm  gives  a  better  estimate  of 
chlorophyll  (and/or  other)  concentrations  than  do  the  present  methods  of  color  index  and  principal  component 
analysis.  The  Optimal  Design  method  has  been  shown  to  be  very  robust  to  changes  in  experimental  conditions  and 
the  models  used  to  determine  Inherent  Optical  Properties.  The  hyperspectral  detector  provides  a  simple  realization 
of  the  OD  algorithm  with  some  sacrifice  in  retrieval  accuracy  because  of  finite  spectr^  width.  A  limitation  of  this 
study,  and  hence  of  the  resulting  optimal  designs,  was  that  the  radiance  spectra  used  as  the  basis  for  the  OD 
computations  were  simulated  and  not  measured.  Another  limitation  of  the  simulated  data  was  that  fluorescence  and 
the  sea  bottom  effects  were  not  considered. 

The  next  step  in  the  solution  of  the  problem  should  be  an  extensive  collection  of  ocean  radiance  spectra  with  a 
hyperspectral  sensor  together  with  field  truth  measurements  of  the  in  situ  optically  active  materials.  This  would 
provide  an  experimental  basis  for  the  comprehensive  statistic  (w(A),  0)  and  allow  the  computation  of  a  really 
optimal  desiga  In  the  future  we  will  apply  the  OD  theory  to  the  problem  of  bottom  imageiy.  Additionally,  we  will 
t^or  the  sensor  system  and  retrieval  aJgorithm  to  match  the  water  of  the  local  area  of  interest  and  optimize  the 
remote  estimation  of  chlorophyll  and  other  optically  active  oceanic  parameters, 
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